Complete next-to-leading-order study on the yield and polarization of T(1S,2S,3S) at 

the Tevatron and LHC 



m 
o 

(N 



Bin Gong, Lu-Ping Wan, Jian-Xiong Wang and Hong-Fci Zhang 
Institute of High Energy Physics, Chinese Academy of Sciences, P.O. Box 918(4), Beijing, 100049, China. 

(Dated: May 6, 2013) 

Based on nonrelativistic QCD factorization scheme, we present the first complete next-to-leading 
order study on the yield and polarization of T(1S, 2S, 3S) hadroproduction. By using the color-octet 
long-distance matrix elements obtained from fits of the experimental measurements on T yield and 
polarization at the Tevatron and LHC, our results can explain the measurements on the yield very 
well, and for the polarizations of Y(1S, 2S, 35), they are in (good, good, bad) agreement with recent 
CMS measurement, but still have some distance from the CDF measurement. 
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Recently, a very important and interesting measure- 
ment on the polarization of T (IS, 25, 35) at the LHC 
was reported by the CMS collaboration [l| with enploy- 
ing improved consideration in the measurement Q while 
that of J/tp has not been worked out yet. It is said 
that measurement on T is easier than that on J/tp at 
the LHC in contrast to the situation at the Tevatron. 
The well-known J/tp polarization puzzle became obvi- 
ous when the CDF measurement at the Tevatron 0] was 
found completely different from the leading-order (LO) 
theoretical prediction in the framework of nonrelativistic 
QCD (NRQCD) which was proposed as a factoriza- 
tion approach on heavy quarkonium decay and produc- 
tion [5[ . Similar jsituation was also found in T production. 



In early works @, [?J , 
at the Tevatron was 
NRQCD prediction [i 



measurement on polarization of T 
aresented, while corresponding LO 
1 does not coincide with it. 



In last six years, there is some very important progress 
in the next-to-leading order (NLO) QCD correction cal- 
culation. The NLO corrections to color-singlet (CS) JJtp 
hadroproduction have been investigated in Refs. [qL flOf. 
its transverse momentum (p t ) distribution is found to be 
enhanced by 2 — 3 order of magnitude at high p t region, 
and its polarization changes from transverse into longi- 
tudinal at NLO [n|. The results are reproduced at LO 
in a new factorization scheme for large pt quarkonium 

production [ll|. The NLO corrections to J/tp produc- 

i [si 3 rsi 

tion via 5-wave color-octet (CO) states (Sq, S[') are 
studied in Ref. [l2| and the corrections to pt distribu- 
tions of both J/tp yield and polarization are small. In 
Refs. [l3|, NLO corrections for Xcj hadroproduction are 
studied. The complete NLO calculation for prompt J/tp 

3 rgi 

hadroproduction (with Pj included) was given by two 
groups [Til . [l5| , and their predictions for p t distributions 
agree with the experimental measurements at the Teva- 
tron and LHC. The calculation for polarization of direct 
J/tp hadroproduction at NLO QCD was presented by two 
groups [3 E3| ■ The complete NLO calculation of the 
polarization for prompt J/tp hadroproduction was com- 
pleted by our group [l8|] last year. It is known that at 
large p t region the logarithm term ln(p 4 / m c ) may ruin fix 
order perturbative expansion and it is resummed in the 



new factorization scheme mentioned above [ill ]. But it 
is unclear that how large the pt region is where fix order 
calculation works well for J/ip case. 

For T hadroproduction, there are studies on the pt dis- 
tribution of yield and polarization for the CS channel at 
QCD NLO % fTql 111 and at the partial next-to-next- 
to-leading order [2(|. NLO QCD correction to p t dis- 
tribution of the yield and polarization for T(1S, 3S) via 
S-wave CO states is presented in Ref. 2l| , and NLO QCD 
correction to pt distribution of the yield for T (15) via all 
the CO states (include 3 Pj) is presented in Ref. [22|. The 
complete NLO study on polarization of T hadroproduc- 
tion has not yet been achieved since there are more com- 
plicated feeddown than charmonium case. However, the 
advantages for study on T are also obvious. Since bot- 
tom is almost three times as heavy as charm, both QCD 

- p 2 ) and v 2 (v is the veloc- 



coupling constant a s (y 4m^ 

ity of heavy quark in the meson rest frame) are smaller 
in bottomonium case. Then perturbative calculation on 
bottomonium should have better convergence in the dou- 
ble expansion of a s and v 2 than that on charmonium. 
Furthermore, it is known that fix order calculation should 
be good enough at intermediate pt region although the 
logarithm term h\(p t /mo) needs to be resummed at large 
Pt region, hence fix order prediction on T hadroproduc- 
tion will be very good for p t up to 60 GcV if that on J/tp 
is very good for p t up to 20 GeV, where 20 GeV is a very 
conservative estimate. In other words, it is expected that 
the theoretical predictions on the polarization and yield 
of T at QCD NLO should be in better agreement with 
experimental measurement than that of charmonium and 
the logarithm term does not need to be resummed for p t 
up to 60 GeV. Therefore, a full study on the polarization 
and yield on T at QCD NLO is a very interesting and 
important task to fix the heavy quarkonium polarization 
puzzle while there are already polarization measurement 
on T (IS, 25, 35) by CMS. In this Letter, we present the 
first complete NLO study on the polarization and yield 
of T(1S, 2S, 3S) based on NRQCD factorization scheme. 

According to the NRQCD factorization formalism, the 
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cross section for hadroproduction of H is expressed as 

da[pp -> H + X] = J2 J d Xl dx 2 G l p Gl 

i,j,n 

x ^ (bb) n X}(<D%), (1) 

where p is either a proton or anti-proton, the indices i, j 
run over all the partonic species and n denotes the color, 
spin and angular momentum states of the intermediate 
bb pair. It can be ^IV^FW 1 and 3 P [ f ] for T, or 3 P^ ] 
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and S[ for XbJ- The short-distance contribution a can 
be calculated pcrturbatively, while the long-distance ma- 
trix elements (LDMEs) (O^) arc fully governed by non- 
perturbative QCD effects. 

The polarization of T is described by three parameters, 
as defined in Ref. 
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where das z st is the spin density matrix of T hadropro- 
duction. In this work, we focus on polarization parameter 
A in helicity frame only. 

To obtain das^S' j similar treatment as in Ref. (l8j 
is taken for both direct and feeddown contributions. 
There are various feeddown contributions in T produc- 
tion, while some of them are ignored in our calculation 
as they are thought to be small. The feeddown contribu- 
tions included in this work are: 

• T(35): no feeddown contribution is included. 

• T(25): feeddown contributions from Y(35) and 
X6j(2P) are included. 

• T(15): feeddown contributions from T(2S, 35) and 
Xbj{lP, 2-P) are included. 

Newly updated Feynman Diagram Calculation pack- 
age [U is used in our calculation. 

In our numerical calculation, the CTEQ6M PDFs [26| 
and corresponding two-loop QCD coupling constant a s 
are used. Branching ratios and masses of involving bot- 
tomonia can be found in Tab. HI Mass of bottom quark 
is set mi, = Mn/2 as an approximation, while Mh is 
the mass of bottomonium H . The CS LDMEs are es- 
timated by using a potential model result [271 ] . which 
gives |i? T (is,2S,3S)(0)| 2 = 6.477, 3.234, 2.474 GeV 3 , and 
|i?^ b(1P2P) (0)| 2 = 1.417, 1.653 GcV 5 , respectively. The 
rcnormalization and factorization scales arc chosen as 
fj, r = fif = rriTi with mr = y (2mb) 2 + p 2 , while the 
NRQCD scale is chosen as = mi,v w 1.5 GcV. The 
ccnter-of-mass energy is 1.8 and 1.96 TeV for Tevatron 
run I and run II, and 7 TeV for the LHC, respectively. 
Various rapidity cuts are chosen according to various ex- 
perimental data, including both central and forward ra- 
pidity region. Besides, a shift pf w pf x (Mh/Mjj 1 ) is 
used while considering the kinematics effect in the feed- 
down from higher excited states. 



In Ref. [22[, feeddown contributions from 3 5[ 8 ' chan- 
nel of XbJ are included in corresponding CO LDMEs of 
T. But when studying the polarization, we have to sep- 
arate Si channels of T and XbJ-, as they have different 
behavior in polarization. Unfortunately, there is still no 
experimental data for XbJ hadroproduction till now, so 
there is no direct clue to determine the CO LDMEs of 
XbJ production. Thus we take them as extra variables in 
our fit, and include experimental data for T polarization 
as well. It is known that the double expansion in a s and 
v is not good enough in the small p t regions. Therefore, 
only data in the region p t > 8 GeV are used in our fit. 

In the fit, we have used the experimental data for _p t 
distribution of the differential cross section by CDF [y], 
LHCb H, CMS II and ATLAS [H, and of the po- 
larization by CDF [3l| and CMS Q. Three fits are per- 
formed for T(35, 25, 15) hadroproduction step by step. 
Here arc some notes: 

• x 2 /d-o.f. = 117/37 is obtained in the fit of T(35). 
In order to express the uncertainty from the three 
CO LDMEs in theoretical predictions correctly, a 
rotation method is performed as Ref. [l~8[ . The cor- 
responding A is (-1.40 ± 1.10, 1.92 ± 0.00, 1.73 ± 
0.00) x 10~ 2 GeV 3 . 

• x 2 /d.o.f. = 88/37 is obtained in the fit of 
T(25) with four CO LDMEs, one of which is 
(0x*°( 2p )( 3 5( 81 )). The rotation is performed with A 
as (24.79 ± 6.23, 1.71 ± 1.83, -2.07 ±0.11, -0.28 ± 
0.01) x 10- 2 GeV 3 , and Xbj{2P) feeddown con- 
tributes a fraction of about 35 ~ 76% as p t in- 
creases in T(25) hadroproduction. 

• In the fit of T(15), we have also included data for 
the fraction of T(15) from XbjO-P) feeddown by 
LHCb [13 and obtained the four CO LDMEs with 
X 2 /d.o.f. = 107/63. Furthermore, the obtained 
Xb(lP, 2-P) fraction in T(15) is consistent with the 
CDF measurement [33[. 

All the fitted CO LDMEs can be found in Tab. HH 
With these LDMEs, our predictions for the differential 
cross section of T hadroproduction are shown in Fig. [TJ 
while those for the polarization are shown in Fig. [2j The 
uncertainty bands in the figures come from errors of the 
LDMEs. 

From the figures we see that the predictions on the 
yield of T hadroproduction can explain the experimental 
data very well with very small uncertainty in a wide range 
of pt at the LHC and Tevatron, while for the polarization, 
things are quite different. For T(35), the production is 

dominated by 5} channel, which results a transverse 
polarization in high p t region and make the theoretical 
predictions far and far away from the experimental data 
as pt increases. And it is obvious that the polarization 
can not be explained at LO in v 2 and NLO in a s if un- 
known feeddown contribution from higher excited bot- 
tomonia is negligible. For T(15, 25), the predictions for 
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FIG. 1: Differential cross section for T hadroproduction at the Tevatron and LHC. From left to right: T(IS'), T(2S) and 
T(35). Rows from top to bottom are corresponding to different experimental conditions of CDF RUN I, CMS, LHCb and 
ATLAS. The experimental data are taken from Refs. [(j [28l - [3C| . 




FIG. 2: Polarization parameter A for T hadroproduction at the Tevatron and LHC. From left to right: Y(IS'), T(2S) and 
T(3S). Rows from top to bottom are corresponding to different experimental conditions of CDF RUN II, CMS (\y\ < 0.6) and 
CMS (0.6 < \y\ < 1.2). The CMS and CDF data are taken from Refs. [J Hp]. 
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9.5 


10.023 
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10.23 
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10.269 



TABLE I: Branching ratios and masses of bottomonia are taken from PDG [25| . 



H 


<o H ( 1 4 81 )) 


<e> H ( 3 s[ 81 )} 


<0"W))/m? 


T(15) 


11.15 ±0.43 


-0.41 ±0.24 


-0.67 ±0.00 


T(2S) 


3.55 ± 2.12 


0.30 ± 0.78 


-0.56 ±0.48 


T(35) 


-1.07 ±1.07 


2.71 ±0.13 


0.39 ±0.23 


Xm,(2P) 




2.76 ±0.67 




X60(1P) 




1.27 ±0.16 





TABLE II: The obtained CO LDMEs for bottomonia produc- 
tion (in unit of 10~ 2 GeV 3 ). 



polarizations can explain CMS data well, but still have 
some distance from the measurement by CDF. From the 
measurement at the LHC, it is easy to see that the p t 
distribution for T(15) is of steepest slope, and that for 
T(2S) is of steeper slope than that for T(35). For the 
CO contribution, we can see that p t distribution of x 5q 
is the steepest one and that of 3 Sf is the most flat one 
while that of 3 Pj is sensitive to the choice of NRQCD 
factorization scale fi\. From the numerical results, we 
obtained that the \bj fceddown contribution in T(15) 
becomes dominant as p t increases, and polarization via 
this channel is slightly transverse polarized, combined 
with the fact that the direct part is dominated by 1 Sq 
channel at small pt region, we find that T(15) is almost 
unpolarized at all the pt range. The situation for T(2S) is 
similar to T(15), but with more \bj feeddown contribu- 
tion at small p t range. Therefore, the \bj feeddown con- 
tribution is very important to explain the experimental 
measurement on polarization. Although the experimen- 
tal measurement on the fraction of \bji)-P) feeddown in 
T(1S) is already used in the fit, it is preliminary with 
large errors. 

It is believed that final physical results are indepen- 
dent of the NRQCD factorization scale n\, but the de- 
pendence does exist when theoretical calculation is trun- 
cated at fix order in the perturbative expansion. And 
this dependence can be found when the detailed arrange- 
ment of NRQCD factorization formula is taken in the 
calculation with P-wave intermediate state involved. So 
a better way to present the final results is to take fi\ de- 
pendence into the consideration of uncertainty. In both 
figures, we have also shown the results with fj,\ = mb and 



pb a = A.qcd- It is found that the /^a dependence is quite 
small for the pt distribution of T yield and polarization 
due to small contribution of T(^pj 8 ') and XbJ\Pj^)- 

In summary, we present the first complete NLO study 
on the polarization and yield of T(15, 25, 35) hadropro- 
duction. Based on the calculation of the polarization and 
yield for both direct and feeddown contributions, eleven 
CO LDMEs are obtained by fitting the experimental data 
at the Tevatron and LHC step by step for T(35, 25, 15). 
With different choices of the NRQCD factorization scale 
//A, we find that fi\ dependence is very small in p t dis- 
tribution of the yield and polarization for T even though 
it could be quite large for J ftp where the P-wave compo- 
nent contributions are very large. For p t distribution of 
T yield, the experimental measurements at the Tevatron 
and LHC can be explained very well in a wide range of 
p t . For T(35), the polarization can not be explained at 
LO in v 2 and NLO in a s if unknown feeddown contri- 
bution from higher excited bottomonia is negligible. For 
T(15, 25), the predictions for polarization can explain 
the CMS data well, but still have some distance from the 
CDF data. 

Further study needs to be considered. The relativis- 
ts corrections to J/ip hadroproduction [H| is negative 
and large in small pt range, and this infers that the rel- 
ativistic corrections to T(35) is the largest one among 
T(15, 25, 35) and detailed study may change the result 
of fit. The uncertainty from the badly known fraction 
of XbJ feeddown in the fits for T(15, 25) could be large 
which is not presented in the plots. With feeddown con- 
tribution of Xfcj(3P), the polarization of T(35) may be 
explained as well. Therefore a further precise measure- 
ment on the fraction of XbJ feeddown or on direct T 
production will be very helpful to fix the polarization 
puzzle. 

We are thankful for help from the Deepcomp7000 
project of the Supercomputing Center, CNIC, CAS and 
also the TH-lA project of NSCC-TJ. This work is sup- 
ported, in part, by the National Natural Science Founda- 
tion of China (No. 10935012, and No. 11005137), DFG 
and NSFC (CRC110), and by CAS under Project No. 
INFO-115-B01. 



[1] S. Chatrchyan et al. (CMS Collaboration) (2012), [2] P. Faccioli, C. Lourenco, and J. Seixas, Phys. Rev. Lett. 
1209.2922. 



5 



105, 061601 (2010), 1005.2601. P. Faccioli, C. Lourenco, 
J. Seixas, and H. K. Wohri, Phys.Rev.Lett. 102, 151802 
(2009), 0902.4462. 

[3] T. Affolder et al. (CDF Collaboration), Phys.Rev.Lett. 

85, 2886 (2000), hep-ex/0004027. 
[4] M. Beneke and I. Z. Rothstein, Phys. Lett. B372, 157 

(1996), hep-ph/9509375. 
[5] G. T. Bodwin, E. Braaten, and G. P. Lepage, Phys. Rev. 

D51, 1125 (1995). 
[6] D. E. Acosta et al. (CDF), Phys. Rev. Lett. 88, 161802 

(2002). 

[7] V. M. Abazov et al. (DO), Phys. Rev. Lett. 101, 182004 

(2008), 0804.2799. 
[8] E. Braaten and J. Lee, Phys. Rev. D63, 071501 (2001), 

hep-ph/0012244. 
[9] J. Campbell, F. Maltoni, and F. Tramontano, Phys. Rev. 

Lett. 98, 252002 (2007). 
[10] B. Gong and J.-X. Wang, Phys. Rev. Lett. 100, 232001 

(2008), 0802.3727. 
[11] Z.-B. Kang, J.-W. Qiu, and G. Sterman, Phys.Rev.Lett. 

108, 102002 (2012), 1109.1520. 
[12] B. Gong, X. Q. Li, and J.-X. Wang, Phys. Lett. B673, 

197 (2009), 0805.4751. 
[13] Y.-Q. Ma, K. Wang, and K.-T. Chao, Phys. Rev. D83, 

111503 (2011), 1002.3987. 
[14] Y.-Q. Ma, K. Wang, and K.-T. Chao, Phys. Rev. Lett. 

106, 042002 (2011), 1009.3655. 

[15] M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 106, 

022003 (2011), 1009.5662. 
[16] M. Butenschoen and B. A. Kniehl, Phys.Rev.Lett. 108, 

172002 (2012), 1201.1872. 
[17] K.-T. Chao, Y.-Q. Ma, H.-S. Shao, K. Wang, and Y.-J. 

Zhang, Phys.Rev.Lett. 108, 242004 (2012), 1201.2675. 
[18] B. Gong, L.-P. Wan, J.-X. Wang, and H.-F. Zhang, Phys. 

Rev. Lett 110, 042002 (2013), 1205.6682. 



[19] B. Gong and J.-X. Wang, Phys. Rev. D78, 074011 

(2008), 0805.2469. 
[20] P. Artoisenet, J. M. Campbell, J. P. Lansberg, F. Mal- 
toni, and F. Tramontano, Phys. Rev. Lett. 101, 152001 

(2008), 0806.3282. 
[21] B. Gong, J.-X. Wang, and H.-F. Zhang, Phys.Rev. D83, 

114021 (2011), 1009.3839. 
[22] K. Wang, Y.-Q. Ma, and K.-T. Chao, Phys.Rev. D85, 

114003 (2012), 1202.6012. 
[23] M. Beneke, M. Kramer, and M. Vanttinen, Phys. Rev. 

D57, 4258 (1998), hep-ph/9709376. 
[24] J.-X. Wang, Nucl. Instrum. Meth. A534, 241 (2004), 

hep-ph/0407058. 
[25] J. Beringer et al. (Particle Data Group), Phys.Rev. D86, 

010001 (2012). 

[26] J. Pumplin, D. Stump, J. Huston, H. Lai, P. M. Nadolsky, 

et al., JHEP 0207, 012 (2002), hep-ph/0201195. 
[27] E. J. Eichten and C. Quigg, Phys.Rev. D52, 1726 (1995), 

hep-ph/9503356. 
[28] R. Aaij et al. (LHCb Collaboration), Eur.Phys.J. C72, 

2025 (2012), 1202.6579. 
[29] V. Khachatryan et al. (CMS Collaboration), Phys.Rev. 

D83, 112004 (2011), 1012.5545. 
[30] G. Aad et al. (ATLAS Collaboration), Phys.Rev D87, 

052004 (2013), 1211.7255. 
[31] T. Aaltonen et al. (CDF Collaboration), Phys.Rev.Lett. 

108, 151802 (2012), 1112.1591. 
[32] R. Aaij et al. (LHCb Collaboration), JHEP 1211, 031 

(2012), 1209.0282. 
[33] T. Affolder et al. (CDF), Phys. Rev. Lett. 84, 2094 

(2000), hep-ex/9910025. 
[34] G.-Z. Xu, Y.-J. Li, K.-Y. Liu, and Y.-J. Zhang, Phys.Rev. 

D86, 094017 (2012), 1203.0207. 



